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y—{ Abstract. We study the geometry of the space of densities Dens(M), which 

is the quotient space Diff(M) /Diff M (M) of the diffeomorphism group of a com- 
pact manifold M by the subgroup of volume-preserving diffemorphisms, en- 
dowed with a right-invariant homogeneous Sobolev i? 1 -metric. We construct 
an explicit isometry from this space to (a subset of) an infinite-dimensional 

^> sphere and show that the associated Euler- Arnold equation is a completely 

integrable system in any space dimension. We also prove that its smooth 
solutions break down in finite time. 

i— i Furthermore, we show that the i^-metric induces the Fisher-Rao (informa- 

tion) metric on the space of probability distributions, and thus its Riemannian 

£^ distance is the spherical version of Hcllinger distance. We compare it to the 

Wasserstein distance in optimal transport which is induced by an L 2 -metric 
on Diff(Af). The H 1 geometry we introduce in this paper can be seen as an 
infinite-dimensional version of the geometric theory of statistical manifolds. 
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1. Introduction 



The geometric approach to hydrodynamics pioneered by V. Arnold j3] is 
based on the observation that the particles of a fluid moving in a compact n- 
dimensional Riemannian manifold M trace out a geodesic curve in the infinite- 
dimensional group Diff^(M) of volume-preserving diffeomorphisms (volumor- 
phisms) of M. Arnold's framework is very general. It includes a variety of non- 
linear partial differential equations of mathematical physics — in abstract form 
often referred to as Euler- Arnold or Euler-Poincare equations. 

With a few exceptions, papers on infinite- dimensional Riemannian geometry, 
including diffeomorphism groups, tended to focus exclusively on either strong 
metrics or weak metrics of L 2 -typej]] The interest in the latter has to do with 
the fact that such metrics often represent kinetic energies as in the case of 
hydrodynamics, see e.g., [UH21E3], or arise naturally in probability and optimal 
transport problems with quadratic cost functions (such as the Wasserstein or 
Kantorovich-Rubinstein distance); see e.g., JSJ EH EH EE] ■ 

On the other hand, in recent years there have appeared a number of interesting 
nonlinear evolution equations described as geodesic equations on diffeomorphism 
groups with respect to weak Riemannian metrics of Sobolev i/ 1 -type, see e.g., 
(U UHl EU] and their references for examples such as the Camassa-Holm, the 
Hunter-Saxton or the Euler-a-equations. 

In this paper we focus on the H 1 metrics both from a differential-geometric 
and a dynamical systems perspective. We will show that they arise naturally 
on (generic) orbits of diffeomorphism groups in the space of all Riemannian 
metrics on M. The main results of this paper concern the geometry of a subclass 
of such metrics, namely, degenerate right-invariant H l Riemannian metrics on 
the full diffeomorphism group Diff(M) and the properties of solutions of the 
associated geodesic equations. The normalized metric is given at the identity 
diffeomorphism by 



It descends to a non-degenerate Riemannian metric on the homogeneous space 
of right cosets (densities) Dens(M) = Diff(Af )/Diff /t (M). Furthermore, it turns 
out that the corresponding geometry is spherical for any compact manifold M. 



Roughly speaking, a Riemannian metric on a Hilbert manifold is strong if it generates a 
complete norm topology on each tangent space, i.e., if the Riemannian norm is equivalent to 
the Hilbert norm; otherwise it is weak. 
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More precisely, we prove that equipped with (1.1) the space Dens(M) is isomet- 
ric to (a subset of) an infinite-dimensional sphere in a Hilbert space. In fact, the 
metric we define on Dens(M) can be viewed as an analog of the metric intro- 
duced by Otto [30J in the theory of mass transport; furthermore the Riemannian 



distance of (1.1), which will be shown to coincide with the (spherical) Hellinger 
distance well-known in probability and mathematical statistics, can be viewed as 
an analogue of the L 2 -Wasserstein distance. Remarkably, it also turns out that 
our metric induces the so-called Fisher-Rao (information) metric and related 
Chentsov-Amari a-connections which have diverse applications in asymptotic 
statistics, information theory and quantum mechanics, see e.g., [HE]. 

We derive the Euler- Arnold equations associated to the general right-invariant 
H 1 metrics which include as special cases the n-dimensional (inviscid) Burgers 
equation, the Camassa-Holm equation, as well as variants of the Euler-a equa- 



tion. In the particular case of the homogeneous if ^metric ( 1.1 ) the Euler- Arnold 
equation has the form 



1.2) p t + u-Vp+y 



12 = ImP 2 ^ 
2P " 2/i(M) 



where u = u(t,x) is a time-dependent vector field on M with divw = p. This 
equation is a natural generalization of the completely integrable one-dimensional 
Hunter-Saxton equation [19] which is also known to yield geodesies on the ho- 
mogeneous space Diff(S' 1 )/Rot(5' 1 ) (the quotient of the diffeomorphism group of 
the circle by the subgroup of rotations), see [20] . 

We prove that the solutions of ( |1.2 ) describe the great circles on a sphere in 



a Hilbert space and in particular the equation is a completely integrable PDE 
for any number n of space variables. The corresponding complete family of con- 
served integrals will be constructed in terms of angular momenta. We point out 
that with a few exceptions (such as the two-dimensional Kadomtsev-Petviashvili, 
Ishimori, and Davey-Stewartson equations) all known integrable evolution equa- 
tions are limited to one space dimension]^] We hope that the geometric properties 



of (1.2) will make it an interesting novel example in this area. 



Furthermore, we show that the maximum existence time for smooth solutions 



of (1.2) is necessarily finite for any initial conditions, with the L°° norm of the 
solution growing without bound as t approaches the critical time. On the other 
hand, the geometry of the problem points to a method of constructing global 



weak solutions of (1.2). We will describe a strategy showing how this can be 
done using a technique of J. Moser. 

It is also of interest to consider the general form of the right-invariant (a-b-c) 
Sobolev H 1 metric on Diff(M) given at the identity by 

(1.3) ((u, v)) := a I (u,v)dp + b / (5u° , 5v ) dp + c / (du b ,dv^) dp, 

Jm Jm Jm 

where u, v G T e Diff (M) are vector fields on M, p is the Riemannian volume 

formj^] b is the isomorphism TM — > T*M defined by the metric on M and a, b 

2 As noted in |14j , derivation and solution of integrable nonlinear evolution partial differen- 
tial equations in three spatial dimensions has been the "holy grail" in the field of integrability 
since the late 1970s. 

3 The volume form /i is denoted by dfx whenever it appears under the integral sign. 
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and c are non-negative real numbers^] We derive the Euler- Arnold equation for 



the metric (1.3) below; while a detailed study of its geometry with the attendant 
curvature calculations will appear in a separate publication [21J. 

Finally, a comment on the functional analytic framework we chose for the 
paper. While our motivations and objectives are directly related to questions in 
analysis and PDE, in order to better present our geometric ideas we will — with 
few exceptions — work primarily with objects (function spaces, diffeomorphism 
groups, etc.) which consist of smooth functions. However, we emphasize that 
when these objects are equipped with a suitably strong topology (for example, 
any Sobolev H s topology with sufficiently large s > n/2 + 1 will do for our 
purposes) then our constructions are rigorously justified in a routine manner. 
We will not belabour this point and instead refer the reader to the papers [T2] 
and [2S] where such questions are considered in greater detail. 

The structure of the paper is as follows. In Section [2] we review the geometric 
background on Euler-Arnold equations and describe the space of densities used 
in optimal transport, as well as reductions of the L 2 - and if x -type metrics on 
Diff(M), its subgroup Diff At (M), and their quotient. 

In Section [3] we introduce the homogeneous if ^metric on the space of densities 
Dens(M) = Diff(M)/Diff At (M) and study its geometry. Generalizing the results 
of [231 [21] for the case of the circle we show that for any n-dimensional manifold 
the space Dens(M) is isometric to a subset of the sphere in L 2 (M,dfi) with the 
induced metric. The corresponding Riemannian distance is shown to be the 
spherical Hellinger distance. 

In Section [4] we describe the relation of the if ^metric to geometric statistics 
and probability. In particular, we show that on the space Dens(M) it plays the 
role of the classical Fisher- Rao metric. In the case M = S 1 we then use it to 
introduce the analogues of dual affine connections generalizing the constructions 
of Chentsov and Amari. 

In Section [5] we study local properties of solutions to the corresponding Euler- 
Arnold equation and demonstrate its complete integrability, as a geodesic flow 
on the sphere. Since for M = S 1 our equation reduces to the Hunter-Saxton 
equation we thus obtain an integrable generalization of the latter to any space 
dimension. 

In Section [6] we turn to global properties of solutions. We derive an explicit 
formula for the Jacobian, prove that solutions necessarily break down in finite 
time and present an approach to construct global weak solutions. 

In Section [7] we derive the Euler-Arnold equation for the general a-b-c metric 



(1.3) and show that several well-known PDE of mathematical physics can be 
obtained as special cases. We also discuss the situations in which some of the 
coefficients a, b, or c are zero. 

Finally, in Section [8] we present a geometric construction which yields right- 



invariant metrics of the type (1.3) as induced metrics on the orbits of the dif- 
feomorphism group from the canonical Riemannian L 2 structures on the spaces 
of Riemannian metrics and volume forms on the underlying manifold M. 



4 We can allow a < as long as ( 1.3 1 is non- negative — e.g., when the first cohomology of M 
is zero. 
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2. Geometric background 

2.1. The Euler- Arnold equations. In this section we describe the general 
setup which is convenient to study geodesies on Lie groups and homogeneous 
spaces equipped with right-invariant metrics. 

Let G be a (possibly infinite-dimensional) Lie group with a group operation 
denoted by (77, £) 1— >■ 77 o £. In our main examples group elements will be diffeo- 
morphisms and the operation will be their composition. We shall use T e G to 
denote the Lie algebra of G, where e is the identity element. For any rj e G the 
group adjoint is the map Ad^ : T e G — > T e G given by the differential 

u — > Ad v u = (L v o R v -i)^ e u 

where L v stands for a left-translation £ 1— )■ rj o £ on the group while i?^ is the 
corresponding right-translation £ i-> £ o 77. The algebra adjoint ad u : T e G — > T e G 
is given by 

d 



a d u 

at 



Ad^( t ), 
t=o 



where ^(t) is any curve in G with 77(0) = e and r)(0) = w. If the group operation 
is composition of diffeomorphisms then in terms of the standard Lie bracket of 
vector fields we have ad M t> = — [u, v]. 

We equip G with a right-invariant (possibly weak) Riemannian metric which 
is determined by an inner product ((■,■)) on the tangent space at the identity 

((R v * e u, R v * e v)) v = ((u, v)) , 

where rj G G and u, v £ T e G. The Euler- Arnold equation on the Lie algebra for 
the corresponding geodesic flow has the form 

(2.1) u t = —B(u,u) = — ad*w, 

with u{0) = uq and where the bilinear operator B on T e G is defined by 

(2.2) ((B(u,v),w)) = ((u,ad v w)). 



R v * e u, 77(0) = e. 



Equation (2.1) describes the evolution in the Lie algebra of the vector u(t) 
obtained by right-translating the velocity along the geodesic 77 in G starting at 
the identity with initial velocity u(0) = u . The geodesic itself can be obtained 
by solving the Cauchy problem for the flow equation 

ch] 
dt 

Remark 2.1. Rewriting Equation ( 2.1[ ) on the dual space T*G in the form 

|(Ad» = 0, ((Ad;u,v)) = ((u,Ad v v)) 

gives a conservation law Ad* , t \u(t) = u expressing the fact that u{t) is confined 
to one and the same coadjoint orbit during the evolution. 



6 



B. KHESIN, J. LENELLS, G. MISIOLEK, AND S. C. PRESTON 



Remark 2.2. Let if be a closed subgroup of G. A right-invariant metric on G 
descends to an invariant (under the right action of G) metric on the homogeneous 
space G/H if and only if the projection of the metric to T^H C T e G, the 
orthogonal complement to H in the group G, is bi-invariant with respect to 
the subgroup H action. (If the metric in G is degenerate along the subgroup 
H then this condition reduces to the metric bi-invariance with respect to the 



if-action, see e.g., [20]. We shall consider the general case in Section 2.4 below.) 
The corresponding Euler-Arnold equation is then defined similarly as long as 
the metric is non-degenerate on the quotient G/H. 

2.2. Examples: equations of fluid mechanics. We list several equations of 
mathematical physics that arise as geodesic flows on diffeomorphism groups. 

2.2.1. Let G = Diff M (M) be the group of volume-preserving diffeomorphisms 
(volumorphisms) of a closed Riemannian manifold M. Consider the right- 
invariant metric on Diff At (M) generated by the L 2 inner product 

(2.3) ((u,v)) L 2 = / (u,v) dp. 

J M 



In this case the Euler-Arnold equation (2.1) is the Euler equation of an ideal 
incompressible fluid in M 

(2.4) ut + V u u = -Vp, divw = 0, 

where u is the velocity field and p is the pressure function, see [3] . In the vorticity 
formulation the 3D Euler equation becomes 

u t + [u, oj] = , where uj = curl w . 

2.2.2. Consider the right-invariant metric on Diff /i (M) given by the H l inner 
product 

((u,v)) H i = / {(u,v) + a 2 {dv 9 , dv")^ dp,. 

JM 

The corresponding Euler-Arnold equation is sometimes called the Euler-a (or 
Lagrangian-averaged) equation and in 3D has the form 

(2.5) uj t + [u, co] = , where u = curl u — a 2 curl Au; 
see e.g. [18]. 



2.2.3. Another source of examples is related to various right-invariant Sobolev 
metrics on the group G = Diff(5' 1 ) of all circle diffeomorphisms, as well as its 
one-dimensional central extension, the Virasoro group. Of particular interest 
are those metrics whose Euler-Arnold equations turn out to be completely inte- 
grable. 

On Diff(S' 1 ) with the metric defined by the L 2 product the Euler-Arnold equa- 



tion (2.1) becomes the (rescaled) inviscid Burgers equation 

(2.6) u t + 3uu x = , 

while the H 1 product yields the Camassa-Holm equation 

(2.7) Ut — Utxx + 3uu x — 2u x u xx — uu xxx = . 
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Similarly, the homogeneous part of the H 1 product gives rise to the Hunter- 
Saxton equation 

(2.8) Utxx + ^ u x u xx + uu xxx = . 

More precisely, in the latter case one considers the quotient Diff(S' 1 )/Rot(5' 1 ) 
whose tangent space at the identity coset [e], i.e., the coset corresponding to 
the identity diffeomorphism, can be identified with periodic functions of zero 
mean. The right-invariant metric on Diff(5 1 )/Rot(5' 1 ) is defined by the H 1 
inner product on such functions 



({u,v))jji = / u x v x dx 
Js 1 

and the corresponding Euler-Arnold equation is given by the Hunter-Saxton 



equation (2.8) 



We also mention that if G is the Virasoro group equipped with the right- 
invariant L 2 metric then the Euler-Arnold equation is the periodic Korteweg-de 
Vries equation 

u t + 3uu x + cu xxx = , 

which is a shallow water approximation and the classical example of an infinite- 
dimensional integrable system. We refer the reader to [20] for more details on 
these constructions. 

Remark 2.3. The Hunter-Saxton equation will be of particular interest to us 
in this paper. In [2U] |2"1] Lenells constructed an explicit isometry between 
Diff(S ,1 )/Rot(S' 1 ) and a subset of the unit sphere in L 2 (S' 1 ) and described the 



corresponding solutions of Equation (2.8) in terms of the geodesic flow on the 



sphere. Although the solutions exist classically only for a finite time they can 
be extended beyond the blowup time as weak solutions, see [25]. In the sec- 
tions below, we shall show that this phenomenon can be established for flows on 
manifolds of arbitrary dimension. 

2.3. The L 2 -optimal transport and Otto's calculus. Given a volume form 
fi on M there is a natural fibration of the diffeomorphism group Diff(M) over 
the space of volume forms of fixed total volume p(M) = 1. More precisely, 
the projection onto the quotient space Diff(M)/Diff M (M) defines a smooth ILH 
principal bundk^] with fibre Diff^(M) and whose base is diffeomorphic to the 
space Dens(M) of normalized smooth positive densities (or, volume forms) 



Dens(M) = \v e 9J\M) : v > 0, / du 

I J M 

see Moser [2H]- Alternatively, let p = dv jd\x denote the Radon-Nikodym deriv- 
ative of v with respect to the reference volume form fi. Then the base (as the 
space of const ant- volume densities) can be regarded as a convex subset of the 
space of smooth functions on M 

M = Ler(M,l >0 ) : f pd/2=l 

I J M 



5 In the Sobolev category Diff(M) -> Difr s (M)/Diff^(A/) is a C° principal bundle for any 
sufficiently large s > n/2 + 1, see [12] . 
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L -optimal transport 



i \ (w-h 1 Wasserstein metric 

Diff,/M) n(il) = Jac (1 (Ti ) 1 



Figure 2.4. T/ie fibration of Diff(M) /i&er Diff M (M) determined by the 
reference density fi together with the L 2 -metrics. The density \i itself corresponds 
to the constant function equal to 1. 



In this case the projection map 7r: Diff(M) — > M. can be written explicitly 
as it{rj) = Jac At (?7 _1 ) where Jac M (i?) denotes the Jacobian of r\ computed with 
respect to /i, that is, i]*fi = Jac A1 (?7)/x. 

The fact that (rj o £)*/i = £*?7*/i implies that 

(2.9) Jac^r? o f ) = (Jac At (ry) o f ) • Jac^(£). 

As a consequence, the projection n satisfies n(ri o £) = 7r(r/) whenever £ G 
Diff M (M), i.e., whenever Jac M (£) = 1. Thus 7r is constant on the left cosets and 
descends to an isomorphism between the quotient space of left cosets to the 
space of densities. 

The group Diff(M) carries a natural L 2 -metric 

(2.10) ((u o rj, v o r])) L 2 = / (u o rj, v o rj) dfi = / (it, i> ) Jac At (?7 _1 ) d// 

where «, i> G T e Diff(M) and 17 G Diff(M). This metric is neither left- nor right- 
invariant, although it becomes right-invariant when restricted to the subgroup 
Diff M (M) of volumorphisms and becomes left-invariant only on the subgroup of 
isometries. Its significance comes from the fact that a curve t — >■ rj(t) in Diff(M) 
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is a geodesic if and only if t — > i](t)(p) is a geodesic in M for each p e 
Following Otto [30J one can then introduce a metric on the base for which the 
projection 7r is a Riemannian submersion: vertical vectors at T^Diff(M) are those 
fields u o 7] with div (pu) = 0, and horizontal fields are of the form V/ o r\ for 
some /: M — » K, since the differential of the projection is tt*(u or]) = — div (pv) 
where p = tt(t])^ 

The induced metric on the base is then 

(2.11) ((a,/3» P = / P (Vf,Vg)dp, 

J M 

where / and g solve div (pV/) = —a and div (pVg) = —f3 with mean- zero 
functions a and /3 considered as elements of the tangent space at p. 



The geodesic equation of the metric (2.10) on Diff(M) is 

dt dt 

i.e., individual particles of the geodesic flow of diffeomorphisms rj(t) move along 
the geodesies in M until they cross (and a smooth solution ceases to exist). In 
Eulerian coordinates, using jWt,s) = u(t,r)(t,x)), the geodesic equation can 
be rewritten as the pressureless Euler (or, inviscid Burgers) equation 

(2.12) d t u + V u u = 

and the induced geodesic equation on the quotient space reads 

(2.13) <9^+I|V0| 2 = O 

where u = V0 and <f> : M — > K is a smooth function]^] Note that, as with any 
Riemannian submersion, if the tangent vector of a geodesic is initially horizontal, 
then it will remain so at later times. In our situation this corresponds to the fact 



that vorticity is conserved so that if a solution of (2.12) is initially a gradient, 
it will always be a gradient (and u = V0 implies that satisfies (2.13)). 

The associated Riemannian distance in Diff(M)/Diff M (M) between two mea- 
sures v and A has an elegant interpretation as the L 2 -cost of transporting one 
density to the other 



(2.14) dist l4 /(z/, A) = inf / dist M (x 1 r](x)) dp 

v Jm 

with the infimum taken over all diffeomorphisms r\ such that rfX = v and where 
distM denotes the Riemannian distance on M; see [S] or [20]. The function dist^ 
is the L 2 -Wasserstein (or Kantorovich-Rubinstein) distance between p and v and 
is of fundamental importance in optimal transport theory. 



Note that if M is flat then (2.10) is a flat metric on Diff(M): a neighborhood of the 
identity is isometric to an open set in the pre-Hilbert space of vector-valued H s functions on 
M equipped with the L 2 inner product. 

7 Otto's construction actually comes from Jac M (ry _1 ), which is important since it is left- 
invariant, not right-invariant. 

8 The projection of the pressureless Euler equation onto the subgroup of volumorphisms 
Diff AI (Af) is the Euler equation of incompressible fluids, so that in a sense fluid mechanics is 
"orthogonal" to optimal transport, as remarked in |30j . 
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2.4. Homogeneous metrics on Diff(M)/Diff M (M). In this section we formu- 
late a condition under which metrics on Diff(M) descend to the homogeneous 
space of densities and describe several examples. The non-invariant L 2 metric 
used in optimal transport, as well as our main example, the right-invariant H 1 



metric (1.1), both descend to the quotient. But other natural candidates, such 
as the right-invariant L 2 metric on Diff(M) or the full H l metric, do not. 

We start with a general observation about right-invariant metrics. Let H be 
a closed subgroup of a group G. 

Proposition 2.5. A right-invariant metric ((■,•)) on G descends to a right- 
invariant metric on the homogeneous space G/H if and only if the inner product 
restricted to T^H (the orthogonal complement ofT e H) is bi-invariant with re- 
spect to the action by the subgroup H, i.e., for any u,v G T^H C T e G and any 
w G T e H one has 

(2.15) ((v , ad^w)) + ((u, &d w v}) = . 

Proof. The proof repeats with minor changes the proof for the case of a metric 
that is degenerate along a subgroup H\ see [20J. In the latter case condition 



(2.15) reduces to bi-invariance with respect to the if-action and there is no need 
to confine to the orthogonal complement T^H. We only observe that in order 
to descend the orthogonal part of the metric must be Ad#-invariant 

({Ad h u, Ad h v)) T ± H = ((u, v)) T ± H = ({u, v}} {h G H) 

for any u, v G T^H. The corresponding condition for the Lie algebra action is 
obtained by differentiation. □ 

Remark 2.6. It can be checked that the condition in Proposition |2.5| is pre- 
cisely what one needs in order for the projection map from G to G/H to be a 
Riemannian submersion, i.e., that the length of every horizontal vector is pre- 
served under the projection. Since the metric on G is assumed right-invariant, 
this condition reduces to one that can be checked in the tangent space T e G at 
the identity. 



Example 2.7. The degenerate right-invariant H l metric (ITT) on Diff(M) de- 



scends to a non-degenerate metric on the quotient Dens(M) = Diff(M) /DifL(M) 
The skew symmetry condition (2.15) in this case will be verified in Section ulbe- 



low (Corollary 7.4). 



Example 2.8. The right-invariant L 2 -metric (2.3) does not verify (2.15) and 
hence does not descend to Dens(M). In fact, if we set u = v = V/ then for any 
vector field w with div w = integration by parts gives 



((ad w v,v)) L 2= / (V v w,v)dfi= / (fVAf,w)d t i = - ((Af)Vf,w)dfi, 

JM JM JM 

where we used the identity div Vv/w = div V^V / — w ■ Af. It is not difficult to 
find / and w such that the above integral is non-zero. For example, we can take 
w to be the divergence-free part of the field (A/)V/ and arrange for a suitable 
/ G C°°(M) so that w ^ 0. 

Similarly, it follows that the full H l metric on Diff(M) obtained by right- 



translating the a-b-c product (1.3) also fails to descend to a metric on Dens(M). 



Note that the c-term in (1.3) does not contribute in this case. 
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Example 2.9. As already pointed out in Section 2.3 the non-invariant L 2 met- 



ric (2.10) descends to Otto's metric on the quotient space whose Riemannian 
distance is the L 2 -Wasserstein distance on Dens(M). This metric is invariant 
under the action of Diff M (M). 

Remark 2.10. One can also consider non-invariant Sobolev H 1 metrics anal- 
ogous to the non-invariant metric ( 2.10[ ) on Diff(M). If the manifold M is flat 



o, 

then (identifying a neighbourhood of the identity in the diffeomorphism group 
with a neighbourhood of zero in a vector space) the energy functional of such a 
metric evaluated on a curve rj(t) will have the form 

(a (d t rj, d t rf) + b {Sd t i], 6d t r)) + c (dd t rj, dd t rj) ) dfidt . 

io J M 

The first term will not be affected by a volume-preserving change of variables. 
However, the terms involving derivatives in the space variable (5 and d) will not 
be conserved in general. This argument can be developed to show that among 
metrics of this type the non-invariant L 2 metric (corresponding to b = c = 0) is 
the only one descending to the homogeneous space of densities Dens(M). 

Table 1. The geometric structures associated with L 2 and H l 
optimal transport. 



Diff(M) 


Diff M (M) 


Dens(M) = Diff(M)/Diff (U (M) 


L 2 -metric 
(non-invariant) 


L 2 -right invariant metric 
(ideal hydrodynamics) 


Wasserstein distance 
(L 2 -optimal transport) 


iP-metric 
(right-invariant) 


Degenerate 
(identically vanishing) 


Spherical Hellinger distance 
(iP-optimal transport) 



3. The ^-spherical geometry of the space of densities 

In this section we study the homogeneous space of densities Dens(M) on a 
closed n-dimensional Riemannian manifold M equipped with the right-invariant 



metric induced by the H 1 inner product (1.1), that is 



(3.1) ((u o rj, v o rffjjji = - I divu-divvdfi 

4 J M 

for any m,d6 T e Diff(M) and rj G Diff(M). It corresponds to the a = c = term 



in the general ia-b-c) Sobolev H 1 metric (1.3) of the Introduction in which, to 
simplify calculations, we set b = 1/4. (We will return to the case of b > in 
Sections [7] and |8}) 

The geometry of this metric on the space of densities turns out to be par- 
ticularly remarkable. Indeed, we prove below that Dens(M) endowed with the 



metric (3.1) is isometric to an open subset of a round sphere in the space of 
square-integrable functions on In the next section, we will show that (3.1) 



9 This construction has an antecedent in the special case of the group of circle diffeomor- 
phisms considered by Lenells [23l [21] • 
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corresponds to the Bhattacharyya coefficient (also called the affinity) in proba- 
bility and statistics and that it gives rise to a spherical variant of the Hellinger 
distance. Thus the right-invariant if 1 -metric provides good alternative notions 
of "distance" and "shortest path" for (smooth) probability measures on M to the 
ones obtained from the L 2 -Wasserstein constructions used in standard optimal 
transport problems. 



3.1. An infinite-dimensional sphere. We begin by constructing an isometry 
between the homogeneous space of densities Dens(M) and an open subset of the 
sphere of radius r 



S™ = \feL\M,dn): / f 2 d^ = r 

I JM 

in the Hilbert space L 2 (M,dfi). 

As before, we let Jac M (?]) denote the Jacobian of r\ with respect to the reference 
form n and let /i(M) stand for the total volume of M. 

Theorem 3.1. The map $ : Diff(M) -»■ L 2 (M,dfi) given by 

$ '. Tj i y f = a/ Jac^r? 

defines an isometry from the space of densities Dens(M) = Diff(M)/Diff At (M ) 
equipped with the H 1 -metric (3.1 ) to an open subset of the sphere C L 2 (M, d/i) 
of radius 



r = y/fi(M) 

with the standard L 2 metric. 

For s > n/2 + 1 the map $ is a diffeomorphism between Diff s (M)/Diff* (M) 
and the convex open subset of (1 H S ~ 1 (M) which consists of strictly positive 
functions on M . 

Proof. First, observe that the Jacobian of any orientation-preserving diffeomor- 
phism is a strictly positive function. Next, using the change of variables formula, 
we find that 

/ $ 2 (r/)(i/i= / J&c^r] dfi = / r]*dpi— dpi — fi(M) 

JM JM JM Jri(M) 

which shows that $ maps diffeomorphisms into S%°. Furthermore, observe that 
since for any £ G Difl^(M) we have 

Jac M (£ o r?)/i = (£ o rj)*n = rj*n = Jac M (ry)//; 

it follows that $ is well-defined as a map from Diff(M)/Diff jU (M). 

Next, suppose that for some diffeomorphisms r\\ and i]2 we have Jac At (r/i) = 
Jac M (?72). Then (771 o r/^ 1 )*/! = /x from which we deduce that $ is injective. 
Moreover, differentiating the formula Jac M (?y)/x = rf \i with respect to rj and 
evaluating at U 6 T^Diff(M), we obtain 

Jac J) » At (?7) = div(U o r/^ 1 ) o r] J&c^r). 
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Therefore, letting it : Diff(M) -»■ Diff(M)/Diff M (M) denote the bundle projec- 
tion, we find that 

((($ o 7r)* TJ (Z7), ($ o n)* v (V))) L 2 = - I (div u o 77) • (div v o 77) Jac^ d/i 

4 Jm 



1 /• 

- / div m • div vdfi = ((U, V)) , 
4 Jm 



4 

for any elements U = uorj and V = v or/ in T 7? Diff(M) where 77 G Diff(M). This 
shows that $ is an isometry. 

When s > n/2 + 1 the above arguments extend routinely to the category of 
Hilbert manifolds modelled on Sobolev H s spaces. The fact that any positive 
function in fl H S ^ 1 (M) belongs to the image of the map $ follows from 
Moser's lemma [29J whose generalization to the Sobolev setting can be found for 
example in [12]. □ 



As an immediate consequence we obtain the following result. 
Corollary 3.2. The space Dens(M) = Diff(M)/Diff |t (M) equipped with the 



right-invariant metric (3.1) has strictly positive constant sectional curvature 
equal to l//i(M). 

Proof. As in finite dimensions, sectional curvature of the sphere equipped 
with the induced metric is constant and equal to 1/r 2 . The computation is 
straightforward using for example the Gauss-Codazzi equations. □ 

It is worth pointing out that the bigger the volume n{M) of the manifold the 
bigger the radius of the sphere and therefore, by the above corollary, the 
smaller the curvature of the corresponding space of densities Dens(M). Thus, 
in the case of a manifold M of infinite volume one would expect the space of 



densities with the if ^metric (3.1) to be "flat." Observe also that rescaling the 



metric (3.1 ) to 



b / div u ■ div v d[i 

' M 



changes the radius of the sphere to r = fi(M). 

3.2. The ^-distance and i^-diameter of Diff(M)/Diff„(M). The right 



invariant metric (3.1) induces a Riemannian distance between densities (mea- 
sures) of fixed total volume on M that is analogous to the Wasserstein distance 
(2.14) induced by the non-invariant L 2 metric used in the standard optimal 



transport. It turns out that the isometry $ constructed in Theorem 3.1 makes 



the computations of distances in Dens(M) with respect to (3.1) simpler than 
one would expect by comparison with the Wasserstein case. 

Consider two (smooth) measures A and v on M of the same total volume 
n{M) which are absolutely continuous with respect to the reference measure [i. 
Let d\/dfi and dv/dfi be the corresponding Radon-Nikodym derivatives of A 
and v with respect to /i. 



14 



B. KHESIN, J. LENELLS, G. MISIOLEK, AND S. C. PRESTON 



Diff(M) 




t Diff(M) 
H l -metric 



H ^optimal transport 



Diff (M) <*>0l) = ^ac^ 



spherical Hellinger 
metric 



Figure 3.3. The fibration of Diff(M) with fiber Diff M (M) determined by the 
reference density fi together with the ^-metric. 



Theorem 3.4. The Riemannian distance defined by the H 1 -metric ( 3.1[ ) between 



measures X and v in the density space Dens(M) = Diff(M)/Diff ft (M) is 



(3.2) dist^A, v) = VMM) arccos ^ \[^ d ^ 

Equivalently, if rj and ( are two diffeomorphisms mapping the volume form \i 
to A and v , respectively, then the H 1 -distance between r\ and ( is 

dist^i (r),() = dist^i (A, v) = a/ /i(M) arccos \ J^j^ J V Jac ^V ' Jac^C df?j . 

Proof. Let f 2 = d\/d[i and g 2 = du/dfi. If A = if /i and v = (*fi then using the 
explicit isometry $ constructed in Theorem 3.1| it is sufficient to compute the 



distance between the functions $(77) = / and $(£) = g considered as points on 
the sphere with the induced metric from L 2 (M, d/i). Since geodesies of this 
metric are the great circles on S%° it follows that the length of the corresponding 
arc joining / and g is given by 

rarccosfl/"/^ 



M 



which is precisely formula (3.2). □ 
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11 



D 



N 



Figure 3.6. A density \i thinly distributed over the manifold and a peaked den- 
sity with the same total volume as \i. 



We can now compute precisely the diameter of the space of densities using 
standard formula 

diam^i Dens(M) := sup {dist^i(A, v) : A, v £ Dens(M)}. 

Proposition 3.5. The diameter of the space Dens(M) equipped with the H l - 
metric (3.1) equals ny/ jj(M)/2. 

Proof. First, observe that for any two densities A and v in Dens(M) we have 




dX dv 
dfj, dfj 



dfj > 0. 



The arc-cosine of this integral is less than 7r/2 which yields ~K\J /j(M)/2 as an 
upper bound for the diameter |*°| 

In order to show that the diameter of Dens(M) is in fact equal to 7Ta/ /i(M)/2 
we construct a sequence of measures such that the distance between them con- 
verges to this limit. Given any large N consider a disk of volume fi(M)/N 
with respect to the reference measure /i. Let be a smooth measure whose 
Radon-Nikodym derivative ffj = du^/d^i is a mollification of iV • xd n , where 
Xd n is the characteristic function of the disk, see Figure 3.6 Note that the total 
volume of un is the same as that of /i. 

We will now estimate the distance between and /i. Since on the disk 
the function f N is approximately equal to y/~N, the integral appearing in the 



argument of arc-cosine in (3.2) can be estimated by 




d[i dv 



N 



dfx d/j 



dfj, 



f%dii 



D 



N 



KM) 
N 



N 



as 



N -> oo. 



10 In particular, this means that the image of DifF(M) under $ is a relatively small subset 
of wh ose p oints can be "seen" from the origin £ L 2 (M,d[i) at an angle less than tt/2, 



see Figure 4.3 
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It follows now that dist^^/i, i/jy) — > 7r-y///(M)/2 as N — > oo which completes 
the proof. □ 

Remark 3.7 (Applications to shape theory). It is tempting to apply the dis- 
tance dist hi to problems of computer vision and shape recognition. 

Given a bounded domain E in the plane (a 2D "shape") one can mollify the 
corresponding characteristic function xe and associate with it (up to a choice 
of the mollifier) a smooth measure Ve normalized to have total volume equal to 



1. One can now use the above formula (3.2) to introduce a notion of "distance" 



between two 2D "shapes" E and F by integrating the product of the correspond- 
ing Radon-Nikodym derivatives with respect to the 2D Lebesgue measure. It is 
not difficult to check that dist^i(z/£;, vp) = tt /2 when ve and vp are mutually 
singular and that dist^^z/^, vp) = whenever they coincide. This works also in 
the case when one of the measures is of delta-type and the other is very thinly 
distributed over a large area. 

In this context it is interesting to compare the spherical metric to other right- 
invariant Sobolev metrics that have been introduced in shape theory. For exam- 
ple, in [22] the authors proposed to study 2D "shapes" using a certain Kahler 
metric on the Virasoro orbits of type Diff(S' 1 )/Rot(S' 1 ), see e.g., [221 EE]- This 
metric is particularly interesting because it is related to the unique complex 
structure on the Virasoro orbits. Furthermore, it has negative sectional curva- 
ture. We refer to [32] for details. 

Example 3.8 (H 1 - gradients on the space of densities). The L 2 -Wasserstein 



metric (2.11) induced on the space of densities was used to study certain dissi- 
pative PDE (such as the heat and porous medium equations) as gradient flow 
equations on Dens(M), see [301 EE]- The analogous computations of the H 1 - 
gradients on Dens(M) simplify due to the isometry with the round L 2 -sphere 
discussed above. 

For instance, let H(p) be a functional of the general form 



H(p) = / h(p) dp 
Jm 

where p G Dens(M) and h is a smooth function on densities. Then one computes 
that grad^ii? = h'(p). 

Here is a quick way to see this. For a small real parameter e and any mean-zero 
function /3 on M we have 

H(p + e/3)= [ h(p + e/3)dp= [ h(p) dp + e [ ti \p)/3 dp + 0{e 2 ) . 
Jm Jm Jm 



By Theorem 3^ one can perform the calculation using the L 2 metric on the 
sphere and identify the variational derivative 5H/ 5p of H with its gradient 
grad^iif so that 

which gives the result. 
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Similarly, for the functional 

F(P) :=\J |Vp| 2 p 



one obtains grad^-iF(p) = — Ap. Observe that in this case the associated gradi- 
ent flow equation 

d t p = -grad^xF 
can be interpreted as the heat equation on densities 

dtp = Ap. 

Thus the Dirichlet functional in the right-invariant H l metric on Dens(M) yields 
the same heat equation as the Boltzmann (relative) entropy functional E(p) = 
J M p log p dp in the L 2 -Wasserstein metric. This provides yet one more relation 
of these two gradient approaches to the heat equation, discussed in [2]. 

Remark 3.9. Alternatively, one can carry out the calculation directly in the 
diffeomorphism group Diff(M) as follows. Let t — > r](t) be a curve in Diff(M) 
with 77(0) = 77 and 7/(0) = v o rj where v = V/ G T e Diff(M) and let p = $(77) 



where $ : Diff(M) — > Dens(M) is the projection, see Figure 3.3 Consider a 
functional H(p) lifted to Diff(M) as an invariant functional. We seek the H l - 
gradient of the functional H(r]) at 77 G Diff(M) in the form Vg o rj for some 
function g : M — > IR. One has 

((giad til H,vo V )) til = — 



t=0 JM 



h(J&cl /2 r](t))dp 

- / //(Jac^ 2 ?]) Jac" 1 / 2 // div v o rj Jac^i] dp . 
2 Jm 



where 

(3.3) t| | t=0 Jac M ?y(t) = (div for/) - Jac M r/. 



After changing variables in the integrand, using (3.1), and projecting with the 
help of $ we obtain 

•^(grad^iif) = ^div (gia.d^H o r/ _1 ) o r] ■ J&cj/ 2 !] = h'(p) 

which coincides with the formula obtained above directly on Dens(M). 

4. Probability and infinite-dimensional geometric statistics 
4.1. Spherical Hellinger distance. The Riemannian distance function dist^i 



on the space of densities Dens(M) introduced in Theorem 3.4 is very closely 
related to the Hellinger distance in probability and statistics. 

Recall that given two probability measures A and v on M that are absolutely 
continuous with respect to a reference probability p the Hellinger distance be- 
tween A and v is defined as 
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As in the case of dist^i one checks that distil (A, v) = y/2 when A and v are 
mutually singular and that dist#(A, v) — when the two measures coincide. It 
can also be expressed by the formula 



dist^(A,z/) = ^2(1- BC{X,u)) 

where BC(X, v) is the so-called Bhattacharyya coefficient (affinity) used to mea- 
sure the "overlap" between statistical samples, see e.g., [9j for more details. 

In order to compare the Hellinger distance dist# e z with the Riemannian dis- 
tance dist^i defined in ( |3.2 ) recall that probability measures A and v are nor- 



malized by the condition A(M) = v{M) = /i(M) = 1. As before, we shall 
consider the square roots of the respective Radon-Nikodym derivatives as points 
on the (unit) sphere in L 2 (M,dfi). One can immediately verify the following 
two corollaries of Theorem 13.11 

Corollary 4.1. The Hellinger distance dist# e z(A, v) between the normalized den- 
sities dX = f 2 dfi and dv = g 2 dfi is equal to the distance in L 2 (M,d/i) between 
the points on the unit sphere f,g G C L 2 (M,dfi). 

Corollary 4.2. The Bhattacharyya coefficient BC(X, v) for two normalized den- 
sities dX = f 2 dfi and dv = g 2 dX is equal to the inner product of the corresponding 
positive functions f and g in L 2 (M,d/j,) 



BC(X,u)= [ Jf-%-dp= [ fgdfi. 
Jm V d ^ d V Jm 

Let < a < 7r/2 denote the angle between / and g viewed as unit vectors in 
L 2 (M,dfi). Then we have 

distil (A, v) = 2sin(a/2) and BC(X,u) = cosa, 

while 

dist^i (A, v) — a — arccos_BC(A, v). 

Thus, we can refer to the Riemannian distance distjji(A,v) on Dens(M) as 
the spherical Hellinger distance between A and v. 

4.2. Fisher-Rao information in infinite-dimensions. It is remarkable that 



the right-invariant H l metric introduced in formulas (1.1) and (3.1) provides 
a convenient framework for an infinite- dimensional Riemannian geometric ap- 
proach to mathematical statistics. Efforts directed toward finding suitable dif- 
ferential geometric approaches to statistics have a long history going back to 
the work of Fisher, Rao [21] and Kolmogorov and continued with Chentsov [5] , 
Efron [13] and Amari-Nagaoka [TJ. 

In the classical approach one considers finite-dimensional families of proba- 
bility distributions on M whose elements are parameterized by subsets E of the 
Euclidean space 

S = {u = v SXt ..., Sk e M : (s l5 . . . , s k ) e E c R k } . 

When equipped with a structure of a smooth fc-dimensional manifold such a 
family is referred to as a statistical model. Rao [31] showed that any S carries 
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dist^i (A, v) 



f = *(A) 



g = <$>{v) 




C L\M,dn) 



Figure 4.3. The Hellinger distance distil (A, v) and the spherical Hellinger 
distance dist#i(A, v) between two points f = $(A) and g = &(v) in . The 
thick arc represents the image o/Diff(M) under the map $. 



a natural structure given by a k x k positive definite matrix 



(4.1) 



M 



d log v d log v 

dsi ds i 



udfj, 



k) 



called the Fisher-Rao (information) metricp] 

In our approach we shall regard a statistical model S as a k-dimensional 
Riemannian submanifold of the infinite- dimensional Riemannian manifold of 
probability densities Dens(M) defined on the underlying n-dimensional compact 



manifold M. The following theorem shows that the Fisher-Rao metric (4.1) is 



(up to a constant multiple) the metric induced on the submanifold S C Dens(M) 



by the (degenerate) right-invariant Sobolev /^-metric (1.1) we introduced orig- 
inally on the full diffeomorphism group Diff(M). 



Theorem 4.4. The right-invariant Sobolev ^-metric (3.1) on the quotient 



space Dens(M) of probability densities on M coincides with the Fisher-Rao met- 
ric on any k-dimensional statistical submanifold o/Dens(M). 



Proof. We carry out the calculations directly in Diff(M). Given any v and w 
in T e Diff(M) consider a two-parameter family of diffeomorphisms (si,^) 
r](si, S2) in Diff(M) starting from the identity r/(0, 0) = e with ^-77(0, 0) = v, 



The significance of this metric for statistics was also pointed out by Chentsov [9]. 
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^77(0, 0) = w and let 



v(si, s 2 ) o r/(si, s 2 ) = (si, s 2 ) and io(si, s 2 ) ° 77(51, s 2 ) = ^^( s i> $2) 

be the corresponding variation vector fields along 7/(i, s). 

If p is the Jacobian of r)(s%, s 2 ) computed with respect to the fixed measure /1 



then (4.1) takes the form 

hw = J ^-^logJac M r](si,s 2 )j — (^log Jac^77(si,s 2 ))jac M r7(si,s 2 )rf/i. 



Recall from (3.3) that 



d 

— Jac M r/(si, s 2 ) = divf (si, s 2 ) o 77(51, s 2 ) ■ Jac M r/(si, s 2 ) 

and similarly for the partial derivative in s 2 . Using these and changing variables 
in the integral we now find 



^Jac^(si, g 2 )^Jac M ry(si, s 2 ) , 
(div v orj) ■ (div w orj) Jac^r] d\i 



M 



div t> • div w dfx = 4((v, u>)) ^1 , 

Af 

from which the theorem follows. □ 



Remark 4.5. Theorem 4.4 suggests that the H 1 counterpart of optimal trans- 
port with its associated spherical Hellinger distance is the infinite-dimensional 
version of geometric statistics sought in [TJ and [H]. 

4.3. AfRne connections and duality on Diff(S' 1 )/Rot(5' 1 ). One of the ques- 
tions posed by Amari and Nagaoka asked for an infinite-dimensional theory of 
so-called dual connections, see e.g. P, Section 8.4. Dual connections arise nat- 
urally in classical affine geometry and have turned out to be particularly useful 
in geometric statistics. 

In this section we describe a family of such connections on the density 
space Dens(M) in the case when M = S 1 which generalize the a-connections of 
Chentsov [9]. We will show that and V*""-* are dual with respect to the 



right-invariant if ^metric (3.1) and derive the (reduced) geodesic equations for 



V^. When a = we recover, as expected, the one-dimensional Hunter-Saxton 



equation (2.8). However, in the case a = — 1 we obtain another completely 
integrable system called the /i-Burgers equation, see [26]. The general case of 
an arbitrary compact manifold M will be described in (21]. 

We refer to [1] for basic facts about dual connections on finite-dimensional 
statistical models and to [20] and [21] for details on the H 1 geometry of the ho- 
mogeneous space of densities Diff(S' 1 )/Rot(S' 1 ). It will be convenient to identify 
the latter with the set of circle diffeomorphisms which fix a prescribed point 

Dens(5 1 ) ~ {77 G Diff(5 1 ) : 77(0) = 0} . 
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Denote A = —d% and given a smooth mean-zero periodic function u define the 
operator A^ 1 by 

nx py pi py 

A~ l u{x) — — I / u{z) dzdy + x / u{z) dzdy. 
Jo Jo Jo Jo 

Let v and w be smooth mean-zero functions on the circle and denote by 

V — v or] and W = w o rj the corresponding vector fields on Dens(5' 1 ). For any 

aGlwe define 



V 



(4.2) 77 -> (V^V)(?7) = (w x v + r< a >(i>, w )) o 
where 

(4.3) r^(v,w) = l -±^A- l d x (v x w x ). 

Proposition 4.6. For each a G R £/ie map a right-invariant torsion- 

free affine connection on Dens(S' 1 ) un£/i Christoffel symbols T^ a \ V*- -* i/ie 
Levi-Civita connection of the H l -metric (3.1), u>/w/e V* -1 -* is flat. 

Proof. The first and last assertions are verified in a routine manner using (14. 21) 



and (4.3). That r^ ^ is precisely the Christoffel symbol of the unique (weak) 
Riemannian H 1 connection can be found e.g. in [21]. □ 

Following [1J we say that two connections V and V* on Dens(S' 1 ) are dual 
with respect to ((•, •)) if 

U{(V, W)) = ((VuV, W)) + ((V, V V W)) 

for any smooth vector fields U, V and W. 

Proposition 4.7. The affine connections and V*- - "-* on Dens(S' 1 ) are dual 
with respect to the H 1 -metric (3.1) for any a6R. 

Proof. It suffices to work with right-invariant vector fields U = uorj, V = v or] 
and W = w or] where u, v and w are periodic mean-zero functions on the circle. 
By right-invariance we have U{{V,W))jji = 0. On the other hand, from (4.2) 
and (3.1) for any a we obtain 

(^Pv,w)) H1 + ((vMu a) w)) H1 = 

= 7/ {vxU+^A~ 1 (u x v x ) x ) x w x dx+ ] - I v x {w x u+ k Y L A~ l (u x w x ) x ) x dx 
4 Jo " 4 Jo 

= \J (^ xM ^ Wx ~~ ^u x v x w x + v x (w x u) x - ^v x u x w x ^dx = 

using integration by parts and the explicit formula for A~ l The general case is 
reduced to the above in a routine way as in e.g. [12]. □ 



From the formula (4.3) we see that the Christoffel symbols r^ a ^ do not lose 
derivatives. In fact, with a little extra work it can be shown that this implies 
that is a smooth connection on the H s Sobolev completion of Dens(S' 1 ) to a 
Sobolev-Hilbert manifold if s > 3/2. Consequently, one establishes the existence 
and uniqueness in H s of local (in time) geodesies of using the methods of 
[28] or ED. 
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We next derive the geodesic equations of (reduced to the tangent space 
at the identity) but will not pursue well-posedness questions here. 

Proposition 4.8. The equation of geodesies of the affine a-connection is 
is 

(4.4) u txx + (2 - a)u x u xx + uu xxx = 0. 

The cases a = and a = — 1 correspond to one-dimensional completely inte- 
grate systems: the HS equation (2.8) and the p-Burgers equation, respectively. 

Proof. The equation for geodesies of on Dens(S' 1 ) reads 

ii + T ( n a \r ] ,r ] )=0 

where is the right-translation of Ti a \ Substituting fj = u o r] gives 

u t + uu x + Ti a) (u,u) =0 



and using (4.3) and differentiating both sides of the equation twice in the x 
variable completes the proof. □ 

Remark 4.9. Dual connections of Amari have not yet been fully explored in 
infinite dimensions. We add here that as in finite dimensions [T] there is a simple 
relation between the curvature tensors of i.e. = (1 — a 2 )R^ where 

R(°) is the curvature of the round metric on Dens(S' 1 ). It follows that the dual 
connections V^ -1 ^ and are flat and in particular there is a chart on Dens(5' 1 ) 
in which the geodesies of the latter are straight lines. This permits to write down 
smooth solutions of the Cauchy problem for (4.4) in the case a — 1, i.e. 

u txx + u x u xx + uu xxx = 0, u(0, x) = u (x) 

(with u (Q) = 0) in the explicit form 

u(t,x)=( e tUoAv) dy) / u 0x (y)e tU0 * iy) dy 



where r/ t (x) = r](t,x) = J* e tu ° x ^dy / e tu ° x ^dy is the flow of the solution u. 
5. The geodesic equation: solutions and integrability 



In the preceding sections we studied the geometry of the i7 1 -metric (3.1 ) on the 
space of densities Dens(M). In this section we shall focus on obtaining explicit 
formulas for solutions of the Cauchy problem for the associated Euler-Arnold 
equation and prove that they necessarily break down in finite time. Furthermore, 
we will show how to view the Euler-Arnold equation as a completely integrable 
Hamiltonian system in any space dimension. 

5.1. Clas sica l solutions of the H 1 Euler-Arnold equation. The formalism 



of Section 2.1 applied to the if ^metric on the group Diff(M) yields the following 
equation. 

Proposition 5.1. The Euler-Arnold equation of the right-invariant H 1 -metric 
(3.1) on the homogeneous space Dens(M) reads 

(5.1) pt + u-Vp + lp 2 = — 2p(M) ' where p = divn. 
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Note that in the special case M = S 1 differentiating Equation (5.1) with re- 
spect to the space variable gives the Hunter-Saxton equation (2.8). We postpone 
the proof of Proposition 5.1 until Section [7] where we derive the Euler- Arnold 
equation in the general case of the right-invariant Sobolev if ^metric (1.3) on 



Diff(M), see Corollary 7.4 



Remark 5.2. The right-hand side of Equation (5.1) is independent of time for 
any initial condition p because the integral J M /rdjj, corresponds to the energy 
(the squared length of the velocity) in the if ^metric on Dens(M) and is constant 
along a geodesic. This invariance will also be verified by a direct computation 
in the proof below. 

Consider an initial condition in the form 

(5.2) p(0,x) = divii (:c). 

We already have an indirect method for solving the initial value problem for 
Equation (5.1) by means of Theorem 3.1 We now proceed to give explicit 
formulas for the corresponding solutions. 

Theorem 5.3. Let p = p(t,x) be the solution of the Cauchy problem (5.1)- 

flow of th 
x. Then 

div uq(x) 



(5.2) and suppose that t !->■ r](t) is the flow of the velocity field u = u(t,x), i.e., 



(5.3) 

where 
(5.4) 



u(t, r](t, x)) where rj(0, x) 

p(t,r](t,x)) = 2k tan ( arctan 



2k 



Kt 



4/i(M) 

Furthermore, the Jacobian of the flow is 



(div uq) 2 dp. 



M 



(5.5) 



Jac M (ry(t, x)) = ( cos Kt + 



divu (x) 
2k 



sin Kt 



Proof. For any smooth real-valued function f(t,x) the chain rule gives 



d 
dt 



df 

(f(t, V (t, x))) = -^(t, r)(t, x)) + (u(t, V (t, x)) , V/(t, r)(t, x)) ) 



Using this we obtain from (|5.1|) an equation for / 
(5.6) 



porj 



dt + 2 ~ J 



-C{t) 



where C{t) = (2/i(M)) 1 j M p 2 dp, as remarked above, is in fact independent of 
time. Indeed, direct verification gives 



p(M) 



dC(t) 
dt 



pp t dp 



M 



div u div u t dp 



M 



= — ( (u, V div u) div u dp — - f (div u) 3 dp = 
JM 2 J M 

where the last cancellation follows from integration by parts. 
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Set C = 2k 2 . Then, for a fixed x G M the solution of the resulting ODE in 



(5.6) with initial condition /(0) has the form 

f(t) = 2k tan (arctau(/(0)/2K) - Kt), 



which is precisely (5.3) 



In order to find an explicit formula for the Jacobian we first compute the time 
derivative of Jac M (??)/! to obtain 

— (Jac^/x) = j t (v*V<) =rf{C VL n) = r]* (div u fj,) = (p o rj) Jac M (r?)/i . 
This gives a differential equation for Jac^t], which we can now solve with the 



help of (5.3) to get the solution in the form of (5.5). □ 



Note that (5.5) completely determines the Jacobian regardless of any "am- 



biguity" in the velocity field u satisfying divw = p in equation (5.1). The 
reason is that the Jacobians can be considered as elements of the quotient space 
Dens(M) = Diff(M)/Diff (U (M). (It will be convenient later to resolve the ambi- 
guity by choosing velocity as the gradient field u = VA _1 p.) 



Remark 5.4 (Great circles on S^°). We emphasize that formula (5.5) for the 
Jacobian Jac^r/ of the flow is best understood in light of the correspondence 
between geodesies in Dens(M) and those on the infinite-dimensional sphere 
established in Theorem |3.1| Indeed, the map 

r ~ ~7 divii (o;) . 
t — > W jaCfj,(r}{t, x)j = cos/tt H smKt 

describes the great circle on the sphere C L 2 (M,dp) passing through the 
point 1 with initial velocity Q* e (uo) = ^divu . 

Having an explicit formula for the Jacobian Jac^r? raises a question whether it 
is possible to construct an associated global flow rj : R x M — > M of the velocity 
field u which could be interpreted as a global (in time) weak solution of equation 



(5.1). We address this question in Section 6.2 



5.2. Integrability of the geodesic flow on a sphere. Next, we examine 
in more detail the Hamiltonian structure of the if^-geodesic flow on Dens(M). 
First, however, we recall the finite dimensional case. 

The geodesic flow on the standard sphere. Consider the unit sphere S 1 ™ -1 C W 1 
given by the equation YTj=\ Qj = 1 with q = (qi, . . . , q n ) G lR n and equipped 
with its standard round metric. The geodesic flow in this metric is defined by 
the Hamiltonian H = YTj=iP 2 j on the cotangent bundle T*S n ~ l . It is a classical 
example of a completely integrable system, which has the property that all of its 
orbits are closed. Indeed, the projections of these orbits from ^S 11 " 1 to S 1 ™ -1 
are the great circles. 



Proposition 5.5. The functions hij = Piqj — Pjq%, l<i<j<non T*IR n (as 
well as their reductions to ^S* 1 ' 1 ) commute with the Hamiltonian H = *YTj=\V 2 j 
and generate the Lie algebra so(n). 
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Proof. The Hamiltonian functions hij in T*IR ?1 generate rotations in the (qi,qj)- 
plane in M n , which are isometries of S n ~ l . These rotations commute with the 
geodesic flow on the sphere and hence {h^, H} = 0. 

A direct computation gives {h^, hjk} = h^, which are the commutation rela- 
tions of so(n). □ 

This example illustrates the so-called noncommutative integrability: the geo- 
desic flow may possess many first integrals (in the case of the sphere there are 
n(n — l)/2 of them), however they need not be in involution. Nevertheless, one 
can proceed as follows, see [7J for details. 

Consider the functions 

Hi := h\ 2) 

H 2 := h\ 2 + h 13 + h\ 3 , 

n n n 2 

# n _i := h\ 2 + • • • + h 2 n _ ln = J2p 2 j q2 j ~ ( Yl q i p i) ' 

j=i j=i j=i 

On the cotangent bundle T*S n ~ l the function H n _i coincides with the Hamil- 
tonian H since YTj=\ Qj = 1 anc l YTj^PiQi = ("the tangent plane equation"). 
We thus have 

Proposition 5.6. The functions {Hi : i = l,...,n— 1} form a complete set 
of independent integrals in involution for the geodesic flow on the round sphere 
S n-i c M n^ that ls ^ Hj y = Qj f or any 1 <i,j <n-l. 

Proof. A routine calculation. □ 

Alternatively, one can consider the chain of subalgebras so (2) C so (3) C ... C 
so(n). Then Hk is one of the Casimir functions for so(k + 1) and it therefore 
commutes with any function on so(A; + 1)*. In particular, it commutes with 
all the preceding functions H m for m < k. They are functionally independent 
because at each step H^ involves new functions h^. 

Remark 5.7. Note that all the Hamiltonians are quadratic in momenta. In 
principle, it should be possible to compute them as limits of the quadratic first 
integrals associated with the geodesic flow on an ellipsoid, as the ellipsoid semi- 
axes cij approach the value 1 with different orders 1 + e % for i = 1, . . . , n; we 
again refer to [7] for more details. 

Remark 5.8. The geodesic flow on an n-dimensional elipsoid is known to be 
a bi-hamiltonian dynamical system and its first integrals can be obtained by 
a procedure similar to the Lenard-Magri scheme, see e.g., [3^1 [27J- On the 
other hand, observe that the limit case of the geodesic flow on the sphere in a 
sense manifests the "highest degree" of integrability: the corresponding Arnold- 
Liouville tori are 1-dimensional and all orbits are closed. 

The geodesic flow on an infinite-dimensional sphere. The finite-dimensional con- 
struction of the first integrals H^ described above allows the following infinite- 
dimensional generalization. 
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Consider the unit sphere C L 2 (M,dp). In fact, choose an orthonormal 
basis {ej}'^ 1 in L 2 (M,dp) (for example, the normalized eigenfunctions of the 
Laplace-Beltrami operator on M) and write 

{oo oo 
3=1 3=1 

The Hamiltonian of the geodesic flow on the sphere is given by the integral 
H = J M p 2 dp where p G L 2 (M, dp). 

As in the previous subsection we define functionals hij = PiQj—Pjqi, 1 < i < j, 
the infinitesimal generators of rotations in the (g«, g^-plane. They now form 
the Lie algebra so(oo) of the group of unitary operators on L 2 and, as before, 
generate an infinite sequence of functionally independent first integrals {Hk}™ =1 
in involution. This sequence corresponds to the infinite chain of embeddings 
So(2) C so (3) C ... C so(oo). 

Alternatively, one can introduce a decreasing sequence of Lie algebras so (oo) D 
so(oo — 1) D so(oo — 2) D ... and construct the corresponding conserved 
quantities as follows. Let p^ 1 ' be the orthogonal projection of p E L 2 (M,dp) 
onto the hyperplane e\. Let p^ be the projection of p onto the hyperplane 
(ei,e2) _L C L 2 (M,dp) of codimension 2, etc. We define the q( k ^s similarly and 
then set ^ 

H (k) = f (p(*))2 d/i /" ( q (k)y d[1 _{ [ q {k) p {k) dj \ 
Jm Jm \J m J 

for k = 0, 1, 2, ... . Note that if ^ coincides with the Hamiltonian if = J M p 2 d/x 
on S^°. Furthermore, we have 

Proposition 5.9. The functionals are in involution, i.e., {H^ k ' , H^ m '} = 
for any k, m = 0, 1, 2, . . . . 

Proof. Observe that the functional can be thought of as a Casimir invariant 
for the action of the unitary group in L 2 (M, dp), while if ^ is a Casimir for the 
action of the unitary group on the subspace (ei, ek) ± C L 2 (M, dp). □ 

Either one of the two sequences of functionals {H k } ( £ =1 or {H^} ! ^ =1 con- 
structed above provides infinitely many conserved quantities for the geodesic 
flow on the unit sphere C L 2 (M,dp). 



Remark 5.10. As we mentioned in Remark |5.8| the geodesic flow on the round 
sphere is bi-hamiltonian and completely integrable in a strong sense. This is 
true in the finite- as well as the infinite-dimensional cases. Nevertheless, it 
may still be of interest to exhibit explicitly the bi-hamiltonian structure of the 



Euler- Arnold equation (5.1) and, in particular, examine how the functionals 
constructed in this subsection correspond to the Lenard-Magri type invariants 
of the Hunter-Saxton equation on the circle M = S 1 . 

6. Global properties of solutions 



Explicit formulas of Theorem 5^ make it possible to give a fairly complete 
picture of the global behavior of solutions to the H 1 Euler-Arnold equation 
on Dens(M) for any manifold M. It turns out for example that any smooth 
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solution of this equation (5.1) has finite lifespan and the blowup mechanism can 
be precisely described. 

On the other hand, the fact that the Jacobian of the associated flow is rea- 
sonably well-behaved suggests a way of constructing global weak solutions. By 



the result of Moser [29] the function on the right side of formula (5.5) will be 
the Jacobian of some diffeomorphism as long as it is nowhere zero. In the first 
two subsections below we shall explore the possibility of constructing a smooth 
flow whose Jacobian may vanish at some points (as inevitably happens to the 
solutions past the blowup time). In Section 6.3 we discuss Fredholm properties 
of the associated Riemannian exponential map on Dens(M). 

6.1. Lifts to the diffeomorphism group. First, we note that there can be no 



global smooth (classical) solutions of the Euler- Arnold equation (5.1). As in the 
case of the one-dimensional Hunter-Saxton equation all solutions break down in 
finite time. 

Proposition 6.1. The maximal existence time of a (smooth) solution of the 
Cauchy problem (5.1)-(5.2) constructed in Theorem 5.3 is 



(6.1) 



< T n 



71 1 

1 — arctan I — inf div un(x' 

2k k \2k x£M 



Furthermore, as t /* T max we have \\u(t)\\ci /* oo. 



Proof. This follows at once from formula (5.3) using the fact that divu = p. 
Alternatively, from formula (5.5) we observe that the flow of u(t, x) ceases to be 
a diffeomorphism at t — T max . □ 

Observe that before a solution reaches the blow-up time it is always possible 
to lift the corresponding geodesic to a smooth flow of diffeomorphisms using the 
classical construction of Moser [29J: 

Proposition 6.2. There exists a family of smooth diffeomorphisms rj(t) in 
Diff(M) satisfying (5.5), i.e. such that Jac /i (r^(t)) = ip(t) where 



(6.2) <p{t,x) 
provided that < t < T max . 



divu (x) . \2 

cos nt H sin Kt 

2k J 



Proof. Integrating <p(t,x) over the manifold gives 

sin 2Kt 



(p(t, X) djj = yU(M) COS Kt + 



M 



2k 



div Uq(x) djj, 



+ 



M 

sin 2 Kt 



Ak? 



(div u Q (x)) 2 dfi = n(M) 



M 



where we used formula (|5.4|) for k and the fact that the integral of the divergence 

d(p 



of Mq vanishes. It follows that 



M 



dt 



(t, x) dp = 



2* 
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which allows one to solve the equation Af = —dip/dt for any fixed time t and 
thus produce a function of two variables, e.g., 

(6.3) f(t,x) — fi(x) cos 2nt + f2(x) sin 2nt, 

where fx and /2 are smooth functions on M. 

Since for any t in [0, T max ) the function <p(t, x) is strictly positive one can define 
a time-dependent vector field by X(t,x) = V f(t,x)/ip(t,x). Letting t >->■ £(t) 
denote the flow of X starting at the identity and using Cartan's formula and 
the definition of / we now compute 

since Cx{<£n) = div (ipX)fi. Observing that ip(0) = 1 and £(0) = e we conclude 
that = n for any < t < T max . Since each £(t) is a diffeomorphism, 

letting rj(t) denote its inverse we find that rffi = ip^i, from which it follows that 
Jac M (^(t,x)) = <f(t,x) as desired. □ 



The method of Proposition 6.2 gives a particular choice of a diffeomorphism 



flow 1] and hence a velocity field appearing in (5.1 ) and satisfying div u = (p. The 
flow must break down at the critical time T max since the vector field X becomes 
singular (when if reaches zero). The difficulty here is that one constructs r\ 
indirectly, by first constructing £ = ?7 _1 , and it is this inversion procedure that 
breaks down at the blowup time T max . 

6.2. Heuristics beyond blow-up time. The geometric picture developed in 
Section [3] helps to gain some insight into the blowup mechanism. 

The only functions in the image of the map $ : Diff (M) — > constructed 



in Theorem 3.1 are those that are positive everywhere and the result of Moser 
[29J on transitivity of the action of diffeomorphisms on densities shows that 
these are all such functions. But everywhere-positive functions form a relatively 
small open subset of and so any great circle that emanates from this set 
will eventually reach the antipodal point corresponding to a function which is 
everywhere negative (and thus not a Jacobian). On the other hand, the Jacobian 
of the flow r) is the square of a function on the sphere (and hence it may still 
make sense). 

Example 6.3. A simple toy model of this phenomenon can be constructed in 
finite-dimensional geometric statistics, see pQ . Consider a sample space of three 
outcomes {a, b, c} with probabilities P(a) = p, P{b) = q and P(c) = r where 
p + q + r = 1. The corresponding statistical model S is two-dimensional (its 
local coordinates are p and q). Using the discrete analogue of the Fisher- Rao 



metric (4.1) the map (p, q) (2^/p, 2y/q, 2y/l — p — q) becomes an isometric 
embedding of S into the sphere S 2 in M 3 of radius 2[jj Consider a geodesic in 
S 2 starting at the uniform distribution p = q = r = | 

(6.4) 7(t) = cost + \/2sint, cost — v^sint, ^ costj . 



12 Probability spaces are mapped into the positive octant of the sphere. 
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Figure 6.4. The picture on the left shows a shaded region on a sphere repre- 
senting the positive octant — the image of the square-root map on the probability 
space. A typical spherical geodesic leaves this region in finite time. The picture 
on the right shows a geodesic in the probability space given by squaring each 
component; when the geodesic hits the boundary it "bounces off the wall. " 



The corresponding probabilities are 

P(a) = \ ( \ cost + v^sint) 2 , P(b) = \ cost - v^sint^ 2 



(6.5) P(c) = I cos 2 1 

Thus "global" geodesies of the Fisher-Rao metric "bounce off the walls" of the 
convex triangle as the individual probabilities approach zero and then increase 



again; see Figure 6.4 



Example 6.5. For the Hunter-Saxton equation on Diff(S' 1 )/Rot(S' 1 ) a similar 
phenomenon was explained in In this case the flow r\ is determined (up to 
rotations of the base point) by its Jacobian. Hence, spherical geometry "selects" 
a unique (weak) global solution r](t) G Diff(S' 1 )/Rot(5' 1 ) which is always smooth. 
If the initial velocity is not constant in any interval, then the singularities of r\ 
are isolated so that rj is a homeomorphism (but not a diffeomorphism past the 
blowup time). This may be interpreted as "elastic collisions" of particles (energy 
is conserved). 

In higher dimensions if the Jacobian is not everywhere positive the situation 
is much more complicated. Nevertheless, in this case it may be possible to apply 
the techniques of Gromov and Eliashberg [16] or Cupini et al. [10J in order to 
construct a map with a prescribed Jacobian. Namely, let <p(t, x) be defined by 
the right side of ( 6.2[ ). Then, for any fixed time t, there is a map r](t, x) of M such 



that Jac^^t, x) = (p(t,x) even if (f(t,x) is zero at some x. Furthermore, using 
for example results of [TU], the map n can be chosen to be C°° smooth provided 



that (p is smooth so that by formula (5.5) one may define a "global flow" in the 



space C°°(M, M)/Diff^(M). It would be interesting to extend Moser's argument 
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to construct a global flow of homeomorphisms out of this flow of maps (past the 
blowup time). 

6.3. The H ^exponential map of Diff(M)/Diff M (M). We next describe the 
structure of singularities of the exponential map of our right-invariant if ^metric 
on the space of densities. 

First, recall from Proposition |3.5| that the diameter of Dens(M) with respect 
to the metric (3.1) is equal to ityj n(M)/2. 



Corollary 6.6. Any geodesic in Dens(M) = Diff(M)/Diff (U (M) starting from 
the reference density is free of conjugate points. 

Proof. Indeed, for any geodesic on a sphere the first conjugate point to the 
initial point is its antipode. However, the great arcs contained in the image of 
the space of densities Dens(M) under the isometry $ constructed in Theorem 



3.1 cannot contain antipodal points along them because the diameter of this 
image is 7ry/I(M)/2. □ 

Using the techniques of [28] one can show that the Riemannian exponential 



map of (3.1) on Dens(M) is a nonlinear Fredholm map. In other words, its 
differential is a bounded Fredholm operator (on suitable Sobolev completions of 
tangent spaces) of index zero for as long as the solution is defined. The fact that 
this is true for the general right-invariant a-b-c metric given at the identity by 



(1.3) on Diff(M) or Diff^(M) also follows from the results of [2B]- More precisely, 



we have the following 

Theorem 6.7. For a sufficiently large Sobolev index s > n/2 + 1 the Rie- 
mannian exponential map of (3.1) on the quotient Diff s (M)/Difffj(M) of the H s 



completions is Fredholm up to the blowup time t = T max given in (6.1) 



The proof of Fredholmness given in [2H] is based on perturbation techniques. 
The basic idea is that the derivative of the exponential map along any geodesic 
t i— )■ r)(t) = exp e (tuo) can be expressed as exp ertno = t^dL^t^it), where W(t) 
is a time dependent operator satisfying the equation 

(6.6) *(*)= / A{ry 1 dr+ [ A(ry 1 B(u , ^(r)) dr 

Jo Jo 

and where A = Ad*Ad,, (as long as t < T max ). If the linear operator w i— > 
B(uq,w) is compact for any sufficiently smooth uq then ^(t) is Fredholm be- 
ing a compact perturbation of the invertible operator defined by the integral 
£ A (r)- 1 dr. In the same way one can check that this is indeed the case for the 



homogeneous space of densities with the right-invariant metric (3.1). We will 
not repeat the argument here and refer to |28j for details. 

Remark 6.8. We emphasize that the perturbation argument described above 
works only for sufficiently short geodesic segments in the space of densities. 
Recall that for the round sphere in a Hilbert space the Riemannian exponential 
map cannot be Fredholm for sufficiently long geodesic because any geodesic 
starting at one point has a conjugate point of infinite order at the antipodal 



point. In the case of the metric (3.1) on the space of densities one checks that 



ll-MO 1 |l / oo as t / T max since it depends on the C 1 norm of r\ via the 
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adjoint representation. Therefore the argument of [28] breaks down here past 



the blowup time as equality (6.6) becomes invalidFl 



7. The general Euler- Arnold equation of the a-b-c metric 
In this section we compute the general Euler-Arnold equation for the a-b-c 



metric (1.3) on the full diffeomorphism group Diff(M). We will be particularly 
interested in the degenerate cases when a = and either b = or c = be- 
cause they lead naturally to nondegenerate metrics on the volumorphism group 
Diff M (M) or the homogeneous space of densities Dens(M) respectively. 

It is convenient to proceed with the derivation of the Euler-Arnold equation 
in the language of differential forms. We begin with a brief comment regarding 
the notation. Recall that the exterior derivative operator d on fc-forms and its 
adjoint 5 are related by 

5(3 = {-l) nk+1 *d* (3, (3eQ k+1 {M), 

where the Hodge star operator * is defined with respect to our fixed Riemannian 
volume element /x on M by the formula a A */3 = (a, (3) \x for any a and (3 in 
fl k (M). As usual, the symbols b and (j = b _1 denote the isomorphisms between 
vector fields and one-forms induced by the Riemannian metric on M. While we 
use d and 5 notations throughout, we will continue to employ the more familiar 
formulas when available. For example, in any dimension we have 5v? = — drvu, 
if n = 1 then dvC = 0, while if n = 3 then du b = * curlw 1 ' and 5dvC = (curl 2 uf 
for any vector field, etc. In particular, in dimension n = 3 the a-b-c metric ( |1.3 ) 
assumes the following special form 

((u, v)) = a I (u,v) dfi + b / div u div v dfi + c / (curl w, curl v) dfi. 
Jm Jm Jm 

In dimension n = 2 we have curlw defined by the formula dvi? = (curlu)/i and 
therefore (curl u, curl v) = curlu curlt> is understood as a product of functions. 



For n = 1 the metric (1.3) simplifies to 



((u, v}} — a uv dx + b u x v x dx . 
's 1 Js 1 



Recall also that the (regular) dual T*DiS(M) of the Lie algebra T e Diff(M) 
admits the orthogonal Hodge decomposition^! 

(7.1) T e *Diff(M) = dn°(M) © 5n 2 (M) © U\ 

where Q h (M) and T-L k denote the spaces of smooth fc-forms and harmonic fc-forms 



on M, respectively; see Figure 7.1 



We now proceed to derive the Euler-Arnold equation of the a-b-c metric (1.3) 
Let A : T e Diff(M) T e *Diff(M) be the self-adjoint elliptic operator 

(7.2) Av = av^ + bd5v b + c5dv b 



13 



It is tempting to interpret this phenomenon as the infinite multiplicity of conjugate points 



on the Hilbert sphere forcing the classical solutions of (5.1 ) to break down before the conjugate 
point is reached. 



14 Orthogonality of the components in (7.1) is established for suitable Sobolev completions 



with respect to the induced metric on differential forms ((a",/3")) 
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dQ°(M) 
6-term 




sn 2 (M) 

c-term 



Figure 7.1. The Hodge decomposition of T*Diff(M) . The b and c terms of the 



metric (1.3) are nondegenerate along the axes indicated in the figure. 



(the inertia operator) so that 



(7.3) 



((u, v )) = / (Au, v) dfi, 

<M 



for any pair of vector fields u and v on M. 



Proposition 7.2. The Euler-Arnold equation of the Sobolev a-b-c metric (1.3) 
on Diff(M) has the form 

Aut = —a ((div u) v 9 + i u du + d(u, w)) — b ((div u) d5u b + di u ddu ) 
(7.4) — c ((div u) ddv 9 + L u d5du b + di u 5du b ) 



where A is given by (7.2). 



Proof. By definition (2.2) of the bilinear operator B, for any vectors u, v and w 
in T e Diff(M) we have 



(7.5) 



((B(u } v),w}} = ((u,ad v w)) = — (Au, [v,w])d[j,. 

'm 



Integrating over M the following identity 

(Au, [v, w)) = (d(Au, w),v) — (d(Au, v),w) — dAu(v, w) 

and using 

/ (d(Au, w), v) dfi — — J (Au, w) div v d/i 

we get 



M 



(u,ad v w)) = J ((div v)Au + d(Au,v) + L v dAu,w) dji. 
'm 
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On the other hand, we have 

((B(u,v),w)) = / (AB(u,v),w) d/i 
Jm 

and, since w is an arbitrary vector field on M, comparing the two expressions 
above, we obtain 

(7.6) B(u, v) = A' 1 ((div v) Au + d(Au, v) + L v dAu) . 



Setting v = u, isolating the coefficients a, b, and c, and using (2.1) yields the 
equation (7.4). The simplification in the b term comes from d 2 = 0. □ 



Remark 7.3. Special cases of the Euler- Arnold equation (7.4) include several 
well-known evolution PDE. 

• For n = 1 and a = 0, we obtain the Hunter-Saxton equation 

• For n = 1 and b = 0, we get the (inviscid) Burgers equation 

Ut + 3uu x = 0. 

• For n = 1 and a = b = 1, we obtain the Camassa-Holm equation 

• For any n when a = 1 and b = c = we get the multi-dimensional 
(right-invariant) Burgers equation 

Ut + V u u + udivu + -V|m| 2 = 0; 

(it is sometimes rewritten using (i u du )^ = V u u — ^V|m| 2 and referred to 
as the template matching equation), see [T7] . 

• For any n and a = 6 = c= lwe get the EPDiff equation 

m t + C u m + m div u = 0, m = v 9 — Av!' , 
see e.g., [IB]. 



Now observe that if a = then the a-6-c metric becomes degenerate and can 
only be viewed as a (weak) Riemannian metric when restricted to a subspace. 
There are three cases to consider. 

(1) a = 0,6 7^ 0,c = 0: the metric is nondegenerate on the homogeneous 
space Dens(M) = Diff(M)/Diff M (M) which can be identified with the 
space of volume forms or densities on M. This is our principal example 
of the paper, studied in Sections [3j [5] and [6j 

(2) a = 0, b = 0, c 7^ 0: the metric is nondegenerate on the group of (exact) 



volumorphisms and the Euler- Arnold equation is (7.10), see Corollary 
17.71 and Remark I7T81 below. 
(3) a = 0, b ytz 0, c 7^ 0: the metric is nondegenerate on the orthogonal 
complement of the harmonic fields. This is neither a subalgebra nor 
the complement of a subalgebra in general and thus the approach of 
taking the quotient modulo a subgroup developed in the other cases 
cannot be applied here. However, in the special case when M is the 
flat torus T n the harmonic fields are the Killing fields which do form 
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a subalgebra (whose subgroup isom(T n ) consists of the isometries) . In 
this case we get a genuine Riemannian metric on the homogeneous space 
Diff(T n )/Isom(T n ). 

In cases (1) and (3) above one needs to make sure that the degenerate (weak 
Riemannian) metric descends to a non-degenerate metric on the quotient. This 



can be verified using the general condition (2.15) in Proposition 2.5 



The corollary below contains Proposition 5T which was stated earlier without 
proof. 



Corollary 7.4. If a = c = then the a-b-c metric (1.3) satisfies condition 
(2.15) and therefore descends to a metric on the space of densities Dens(M). 





The corresponding Euler- Arnold equation is 
(7.7) d div u t + div u d div u + di u (d div u) 

or, in the integrated form, 



(jOj) 

where p = div it. 



Pt + (it, Vp> + \p 2 



2p(M) 



Proof. We verify pTSj ) for G = Diff(M), H = Diff M (M) and &d w v = -[w,v], 
where [•, •] is the Lie bracket of vector fields on M. Given any vector fields u, v 
and w with divw = 0, we have 



{{&d w v, u)) 6 i + {(v, &d w ii)) Hl 



— b / ( div [w, v] div u + div [w, u] div v) dp 
—b / ( ((w, V div v) — (v, V div w)) div u 

+ ((w, Vdivu) — (u, Vdivu>)) divv^j dp 
div w ■ div v ■ div u dp = , 



M 



which shows that Q descends to Diff(M)/Diff M (M)Q 

The Euler- Arnold equation on the quotient can be now obtained from (7.4) 
in the form (7.7). In integrated form it reads 

div u t + w(div u) + \ (div u) 2 = C (t) 

where C(t) may in general depend on time. Integrating this equation over M 
determines the value of C(t). □ 

We now return to the nondegenerate a-b-c metric (a ^ 0) and restrict it to the 
subgroup of volumorphisms (or exact volumorphisms) , cf. Figure 7.1 Observe 
that one obtains the corresponding Euler- Arnold equations with b = directly 
from (7.4) using appropriate Hodge projections. 



Corollary 7.5. The Euler- Arnold equation of the a-b-c metric (1.3) restricted 
to the subgroup Diff A1 (M) has the form 

(7.8) au\ + c5du\ + ai u du^ + CL u d5du^ = dA~ 1 5 (aL u du^ + ci u d5du^) . 



15 



Note that the same proof works for the quotient space G/H = Diff(T")/Isom(T rl ). 
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Proof. Let P denote the orthogonal Hodge projection of (Sobolev completions 
of) the tangent space T e Diff(M) onto the tangent space T e Diff M (M) to the sub- 
group at the identity. In this case the Euler- Arnold equation reads^ 



-PB(u, u) = -PA' 1 (5u b Au + d(Au(u)) + L u dAu) 



where we used (7.6) 



Since 5u b = and since the Hodge projection P = I — dA 1 5 commutes with 
the (inertia) operator A we have 



Au t = —Pi u dAu = —L u dAu + dA l 5t u dAu 



which is an equivalent form of (7.8). 



□ 



One should also mention that (7.8) is also the Euler- Arnold equation when 



the a-b-c metric (1.3) is restricted to the set of exact volumorphisms Diff M ex (M), 
i.e. the subgroup of Diff(M) whose tangent space T e Diff /jex (M ) at the identity 
consists of vector fields v of the form v b = 5(3, where (3 is a 2-form on M. In this 
case the corresponding orthogonal projection is P = 5A _1 rf. 



Remark 7.6. It is interesting to note that the Euler- Arnold equation (7.8) is 
closely related to the H 1 Euler-a equation 



(7.9) 



{v 9 + a 2 8du b ) t + C u {u b + a 2 5du b ) = 0, 



(a e 



which was proposed as a model for large-scale motions by Holm, Marsden and 
Ratiu [TBI. 



There is also a "degenerate analogue" of the latter equation which corresponds 
to the case where a = b = 0: 



Corollary 7.7. The Euler- Arnold equation of the right-invariant metric (1.3) 
with a = b = on the subgroup of exact volumorphisms is 

(7.10) 



5du\ + PC u (6du b ) = 0, 



where P is the orthogonal Hodge projection onto 5Q 2 (M). 
Proof. Set a = b = in Equation (7.4) and project onto 5Q 2 (M). 



□ 



Remark 7.8. Observe that the nonlinear PDE in (7.10) has a simpler structure 



than the Euler-a equation (7.9). In fact, one can view equation (7.10) as a formal 



limit of (7.9) when a /* oo. Furthermore, it seems likely that establishing global 



well-posedness for either equation requires overcoming very similar technical 
obstacles and we expect that standard techniques will be sufficient to obtain 



global well-posedness of (7.10) in the two-dimensional case. 



16 Similar computations in the general case of a projection of an Euler-Arnold equation to 
a subgroup can be found in [28] or [37] , 
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8. The space of metrics and the diffeomorphism group 

Our aim in this section is primarily motivational. Apart from the fact that 
the Euler- Arnold equations of H 1 metrics yield a number of interesting evolution 
equations of mathematical physics discussed above there is also a purely geomet- 
ric reason to study them. Below, we show that right-invariant Sobolev metrics 
of the type studied in this paper arise naturally on orbits of the diffeomorphism 
group acting on the space of all Riemannian metrics and volume forms on M. 
First we collect a few basic facts in Section 18.11 



8.1. Weak L 2 structures on the spaces of metrics and volume forms. 

Our main references for the constructions recalled in this subsection are [HI [HJ 

Given a compact manifold M consider the set Met(M) of all (smooth) Rie- 
mannian metrics on M. This set acquires in a natural way the structure of a 
smooth Hilbert manifoldJHl 

The group Diff(M) acts on Met(M) by pull-back g \— y V g (i]) = rfg and there 
is a natural geometry on Met(M) which is invariant under this action. If g is a 
Riemannian metric and A, B are smooth sections of the tensor bundle S 2 T*M, 
then the expression 



U) ((A,B)) g = [ Tr g (AB)dp g 

J M 



defines a (weak Riemannian) L 2 -metric on Met(M). Here \x g is the volume form 
of g and the trace is defined by Ti g (AB) : = Tr (g~ 1 Ag~ 1 B) . This metric is 
invariant under the action of Diff(M), see [TT] . 

The space Vol(M) of all (smooth) volume forms on M also carries a natural 
(weak Riemannian) L 2 -metric 

4 f da dB , 
8.2 ((a,P)) u = - / -r-rdv, 

n J M dv dv 

where v G Vol(M) and a, B are smooth n- forms and which appeared already in 
the paper [I5]rj It is also invariant under the action of Diff(M) by pull-back 

There is a map S: Met(M) — y Vol(M) which assigns to a Riemannian metric 
g the volume form \x g . Its derivative S* 9 : T 9 Met(M) — y T=( s )Vol(M) is 

Z* g (A) = Ti s(g) A = Tr(g- l A). 



One checks that H is a Riemannian submersion in the normalization of (8.2). 
Furthermore, for any g in Met(M) there is a map i g : Vol(M) x {g} — y Met(M) 
given by 



which is an isometric embedding. 



17 Indeed, the closure of C°° metrics in any Sobolev H s norm with s > n/2 is an open 
subset of H S {S 2 T*M). 

18 The space Vol(M) of volume forms on M contains the codimension 1 submanifold 
Dens(M) C Vol(M) of those forms whose total volume is normalized. 
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For any fi G Vol(M) the inverse image Met /Lt (M) = can be given a 

structure of a submanifold in the space of Riemannian metrics whose volume 
form is \i. Its tangent space at g consists of symmetric bilinear forms A with 
Tig(A) = 0. The metric (8.1) induces a metric on Met At (M), which turns it into 
a globally symmetric space. 

The natural action on Met M (M) is again given by pull-back by elements of 
the group Diff M (M). The triple of Met(M) (total space) together with Met^(M) 
(fibre) and Vol(M) (base) forms a trivial fiber bundle, even though (8.1) is not 
a product metric. 

Sectional curvature of the metric (8.1) on Met(M) was computed in [15] and 
found to be nonpositive. The corresponding sectional curvature of Met At (M) is 
also nonpositive. On the other hand, the space Vol(M) equipped with L 2 -metric 
(8.2) turns out to be flat. 

We now explain how these structures relate to our paper. 



8.2. Induced H 1 metrics on orbits of the diffeomorphism group. Ob- 
serve that the pull-back actions of Diff(M) on Met(M) and Vol(M) (and simi- 
larly, the action of Diff M (M) on Met M (M)) leave the corresponding metrics ( |8.1 ) 



and (8.2) invariant. This allows one to construct geometrically natural right- 
invariant metrics on the orbits of a (suitably chosen) metric or volume form. 

We first consider the action of the full diffeomorphism group Diff(M) on the 
space of Riemannian metrics Met(M). 

Proposition 8.1. If g G Met(M) has no nontrivial isometries then the map 



V g : Diff(M) — > Met(M) is an immersion and the metric (8.1) induces a right- 
invariant metric on Diff(M) given at the identity by 

(8.3) ((u,v)) = ({£ v g,£ v g)) g 

for any T e Diff(M). 

Proof. First, observe that the differential of the pull-back map V g {rj) with respect 
to 7] is given by the formula 

*Pg*r,( vo v) = V*{£vg), 

for any v G T e Diff(M) and i] G Diff(M), where C v stands for the Lie derivative. 
If g has no nontrivial isometries then it has no Killing fields and therefore the 
differential V g * is a one-to-one map. 

It follows that the induced Riemannian metric on the orbit through g can be 
computed at any r\ G Diff(M) as 

((u o V ,vo v )) v = ({V g *{u o r)),V g *(v o r)))) v * g 
= ((r]*C u g,7]*C v g)) rg 

for any u, v G T e Diff(M). The induced metric is necessarily right-invariant. In 
fact, writing out explicitly the pull-back in the form rfh = Drj (h o rf)Drf for 
any h in Met(M) and using the definition of Tr 5 and the change of variables 
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formula one readily checks that 

((r]*C u g,r]*C v g)) rg = / Tr v * g (rfC u grfC v g) d/i rg 

JM 

= J Tr [{rf g)" 1 ^* C u g {rf g)" 1 ^* C v g^j \/det(r}*g) dx 



Tig (C u gC v g)dn 9 

M 

= ((C u g,C v g)) g 

holds for any t] e Diff(M), which implies right-invariance. □ 

Remark 8.2. More generally, if g has non-trivial isometries, then the above pro- 
cedure yields a right-invariant metric on the homogeneous space Diff(M) /Iso 5 (M); 
see the diagram (8.7) below. 



Remark 8.3. The inner product (8.3) can be rewritten in terms of the exterior 
differential d and its adjoint 5 so that at the identity we get 

((u, v)) = / (C u g, C v g) dfi 



U) 



M 

(du^dv^dfi + A / (5u\5v b )diJ,~A / Ric( u, v) dfx, 

M JM JM 



for any vector fields u, v e T e Diff(M) and where Ric stands for the Ricci curva- 
ture of M. If the metric g is Einstein then Ric(w, v) = X(u, v) for some constant 
A and the induced metric in (8.3 ) becomes a special case of the a-b-c metric ( 1.3 ) 
with a = —AX, 6 = 4 and c = 2. 

In exactly the same manner we obtain an immersion of the volumorphism 
group Diff„(M) into Met /i (M). 

Proposition 8.4. If g e Met M (M) has no nontrivial isometries then the map 
V g : Diff /i (M) — > Met M (M) is an immersion and (8.1) restricts to a right-invariant 
metric on Diff^(M). 



In this case (8.4) yields the formula 

(8.5) ((u,v)) = 2 (du^dv^dn-A Ric{u,v)dfi, 

JM JM 

for any u, v e T e Diff At (M). 

Finally, we perform an analogous construction for the action of Diff(M) on the 
space of volume forms Vol(M). In this case the isotropy subgroup is Diff^(M) 
and we obtain a metric on the quotient space Diff(M)/Diff At (M). 

Proposition 8.5. If fi is a volume form on M then the map \ Diff(M) — > 
Vol(M) defines an immersion of the homogeneous space Dens(M) into Vol(M) 
and the right-invariant metric induced by (8.2) has the form 

4 f 

(8.6) \{u,v)) = ft = - divu- divvd/i. 

n J M 
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Proof. The differential of the pullback map is 

for any v G T e Diff(M) and rj G Diff(M). Using formula (8.2) and changing 
variables as before we obtain 



n 



n 



(divw o rj) ■ (divf o rj) d(rj*n) 
div u ■ div vdfi = ((u, v)) , 



M 



so that the induced metric is again right-invariant. 



□ 



The three immersions described in Propositions 8.1 , 8.4 and 8J3 can be sum- 
marized in the following diagram. 



Iso s (M) 



(8-7) Iso 9 (M) 



cmb 



-> Diff M (M) Diff M (M)/Iso 3 (M) — ^ Met M (M) 



emb 



emb 



-)■ Diff(M) Diff(M)/Iso 5 (M) 



proj 



cmb 



Met(M) 



Diff M (M) Diff(M) 



-> Diff(M)/Diff„(M) Vol(M) 



The first three terms of each row in (8.7) form smooth fiber bundles in the ob 



vious way. The third column is a smooth fiber bundle since Iso s (M) C Diff M (M). 
The fourth column is a trivial fiber bundle which already appeared in [15] . 

Remark 8.6. While curvatures of the spaces Met(M), Met M (Af) and Vol(M) 
have relatively simple expressions, the induced metrics above on the correspond- 
ing homogeneous spaces 

Diff(M)/Iso s (M), Diff M (M)/Iso g (M) and Diff(M)/Diff M (M) 

turn out to have complicated geometries with one notable exception: the im- 
mersion of Diff(M)/Diff M (M) into Vol(M). This immersion has codimension 
one and, as we discussed throughout this paper, maps the whole of Diff(M) 
into the space Dens(M) C Vol(M) of volume forms which integrate over M to 
the same constant volume as \i. On the other hand, the second fundamental 
form of the immersion of the quotient space Diff(M)/Iso 5 (M) into Met(M) is 
complicated (in this case codimension is infinite). One can show that the sec- 
tional curvature of Diff(M)/Iso g (M) in the induced metric assumes both signs, 
see (ST]. 
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